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Asher Yahalom 
The Racah Institue of Physics, Jerusalem 91904, Israel 



in 

(N 
(N 



> 
O 
oo 
O 

o 
in 

O 



1 



Abstract 

I find conditions under which the "Weak Energy Principle" of Katz, Inagaki and Yahalom 
(1993) gives necessary and sufficient conditions. My conclusion is that, necessary and 
sufficient conditions of stability are obtained when we have only two mode coupling in the 
gyroscopic terms of the perturbed Lagrangian. To illustrate the power of this new energy 
principle, I have calculated the stability limits of two dimensional configurations such as 
ordinary Maclaurin disk, an infinite self gravitating rotating sheet, and a two dimensional 
Rayleigh flow which has well known sufficient conditions of stability. All perturbations 
considered are in the same plane as the configurations. The limits of stability are identical 
with those given by a dynamical analysis when available, and with the results of the 
strong energy principle analysis when given. Thus although the "Weak Energy" method 
is mathematically more simple than the "Strong Energy" method of Katz, Inagaki and 
Yahalom )1993) since it does not involve solving second order partial differential equations, 
it is by no means less effective. 

Key words: Energy variational principle; Self-gravitating systems; Stability of fluids. 
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I. Introduction 



The main purpose of this work is to study the stabihty features of weak variational 
principle of barotropic fluid dynamics (Katz, Inagaki and Yahalom (1993) from now on 
paper I). This principle has the advantage that all equations of barotropic fluid dynamics 
are derived from one Lagrangian, i.e. Euler equations and continuity equation. And 
therefore stability analysis in this formalism does not involve solving partial differential 
equations like in the dynamical perturbation method or the strong variational principle 
(see paper I). However, it has the disadvantage of lacking standard form. This Lagrangian 
contain only a term depending on the degrees of freedom (a potential) and a term linear in 
time derivatives (gyroscopic term), but lacks a term quadratic in time derivatives (kinetic 
term). The above fact is the cause of some novel stability features I discuss below. 

The plan of this paper is as follows, section II contains a quick introduction to the 
weak variational principle of barotropic flows given in (paper I), some modiflcations are 
introduced for both 3-D and 2-D cases . In section III I overview the stability theory of a 
system which does not contain terms quadratic in time derivatives (gyroscopic system), a 
comparison between stability predictions using the strong variational principle and weak 
variational principle is also given. 

In section IV I introduce the general form of the second variation of a barotropic flow 
potentail. Later a few illustrations of stability analysis using the weak energy principle 
in two dimensional flows is given. In section V I give a general formalism of uniformly 
rotating galactic flows. In section VI I study a speciflc case of the above flows, the uniform 
rotating sheet which was flrst analyzed by Binney J. & Tremaine S. (1987). In section VII I 
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analyse the stability of another galactic model known as the "Maclaurin disks" which also 
appears at Binney J. & Tremaine S. (1987). This model is analysed dynamicaly in (Binney 
J. & Tremaine S. (1987)) by solving three equations (Euler equations in two dimensions 
and the continuity equation). It is also analysed in (Yahalom, Katz & Inagaki 1994 from 
now on paper II) using the strong energy method and solving one equation (continuity 
equation). Here it is analysed without solving any equation, but merely diaganolizing the 
appropriate potential. In section VIII I analyse the two dimensional Rayleigh flows , where 
sufficient conditions of stabilty were found by Lord Rayleigh (1880). 

II. The Weak Variational Principle 

The weak variational principle is derived in paper I from the Lagrangian (5.16) com- 
bined with (6.12) : 



v^ + e + l^ ]]pd^x + b-{P-Po) + na- (J - Jo) {11. 1) 



We use the following notations: for the positions of fluid elements r or (x^) — {x,y,z)] 

— * 

iC, L, ... = 1, 2, 3 the density of matter is p; v is the velocity fleld in inertial coordinates, h 
and Vic serve as Lagrange multipliers with respect to linear and angular momenta respec- 
tively given by: 

vpd^x {II.2.a) 



and 

J= J rx vpd^x (I/.2.6) 

e{p) is the speciflc internal energy of the barotropic fluid, related to the pressure and the 
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specific enthalpy: 



P dp 
$ is the internal gravitational potential given by: 



P 8f 
e{p) = h-- P = P^^ {11-^) 



^ = -G I d-'x'. (II A) 



In our formalism v has a Clebsch form: 



And p is given by: 



V = aV/S + Vu. (//.5) 



P=^7 r- [II-Q) 



d{x,y,z) 



w is defined by: 

a + w-Va = $ + w-Wp = fi + w ■ V p, = 0. (//.7) 

Hence our system is described by four trial functions a, f3, p, u and two Lagrange multi- 
pliers b and ^Ic- Taking the variation of Lagrangian (II. 1) with respect to the four trial 
functions gives the Euler and mass conservation equations in moving coordinates. Taking 
the variation of Lagrangian (II. 1) with respect to the two Lagrange multipliers gives the 
fixation of linear and angular momentum. Taking the variation of the potential part of 
(II. 1) gives the equation of a stationary barotropic fluid. Further details can be found in 
paper I. I will know modify the Lagrangian slightly. Rewriting the first part of equation 
(II. 1) using equation (II. 5) : 



/ 



w ■ vpd^x = y ^ • («^/^ + ^i^)pd^x = J {-aPp - vV ■ {pw))d^x (//.8) 



5 



where in the second equahty we integrated by parts and neglected boundary terms (we 
also assumed that v is single valued) and used equation (11.7). Thus we obtain: 



j w ■ vpd^x — J 



w ■ vpd x — {—aPp + i'p)d x= {—aP — i')pd x 



{11.9) 



where in the second equality we added a full time derivative. Inserting equation (II. 9) into 
equation (II. 1) and using equation (II. 5) we obtain: 



L = 



-a$-u- -{aV(3 + Vu)^ + e + -$ 

Zi Zi 



pd^x+b-{p-Po)+nc-{j-Jo). {11.10) 



In 2-D flows the formalism is modified slightly. Euler equations have only two components, 
linear momentum has two components P. Angular momenntum J and vorticity u = 
V X V • 1 ^ have one component. The density p is given by: 



P = J:5d{z) = X{a)^^6D{z). 



(//.II) 



d{x,y) 

5d is diracs delta, S is the surface density, and A is a function of a depending on the flow 
under consideration. The internal energy e and pressure P are considred to be functions 
of E. The Lagrangian of a 2-D flow does not depend on fj, and is given by: 



1 ^ ^ 1 ' 

-a$ -V- -{aV/3 + Vi^)^ + £ + 

Zi z 



j:d^x + b-{p-Po) + nc{J-Jo)- {11-^2) 



For Further details see paper II. 

III. Stability of Gyroscopic Systems 

It is our purpose now to derive from the functional form of the potential term in the 
Lagrangian (11.10) the conditions under which a given stationary conflguration becomes 
stable or unstable. To do this we study a system with N degrees of freedom that is described 
by the Lagrangian: 

L = G-V = bi{q)q'-V{q) {III. I) 
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(a summation agreement is assumed) the energy will become 

E^V{q). (111.2) 

This means that the system can propogate only on a manifold of equipotential, if the 
equipotential manifold is zero dimensional as in the case of a maximum or a minimum of 
the potential than the system is bound to remain on this point i.e. we have a suflBcient 
condition of stability. This can also be deduced from the following argument, suppose 
we distort our equilibrium configuration given hy — qQ,qQ — and f^|go = slightly 
such that initial conditions and energy are slightly diflFerent. Since this is done around 
equilibrium there is no first order contribution to the variation of energy and we obtain: 

= -Q^^UMMj = VijSqiSqj. (111.3) 

From equation (III. 3) it is easy to see that if the matrix Vij has a definite sign (i.e. we 
have a minimum or a maximum) , the system is bounded (i.e. stable). Let us now take 
the dynamical point of view, writing down the linearized Euler-Lagrange equations of this 
system we obtain: 

VijSqi = (|^ - ^)\qoSqi = bijSqi. (III.A.a) 
The above expression can be derived from the perturbed Lagrangian: 

5^L = 5^G - S^V = hjSqiSqj - VijSqiSqj (III.A.b) 

Now suppose that: 

5qi oc e'^* (Ill.b) 
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this means that: 

VijSqi = LjbijSqi. (///.6) 
in order to obtain the eigen frequencies we must solve the equation. 

det\Vij - ubijl = 0. {1 1 1 -7) 



Assume without the loss of generality that Vij is diaganolized. Taking a two dimensional 
perturbation (a two mode coupling in 5^G) we have: 



2 ^11^22 /TTTQ\ 

u = -2 — (///.8) 

in this case a definite sign is a necessary and sufficient condition of stability as can be 
clearly seen from equation (III. 8). Moving on to the third dimensional case: 

2_ ^11^22^33 (TTTa\ 

here a definite sign is not necessary, take for example: 

Vii < 1/22 < ^33 > with 6^2^33 + hl^V22 + &23^11 > 0- (///.lO) 

Thus we conclude that a neseccary and sufficient condition of stability appear only in the 
case where the gyroscopic term couple only two modes, for more complex coupling the 
energy criterion is only sufficient. 

///. 1 Strong systems derived from Gyroscopic systems 

In this section we diverge from the main topic of this paper in order to study the 
relations between weak and strong variational principles (paper I). Observing closely the 



Lagrangian (11.10) we see that it contains two kinds of variables, i.e. we have N + M 
variables Qiji^a (we take i,j,k e [l..iV] and a,b,c e [1..M]) , the Lagrangian has the 
following structure: 

The equations of motion of this system are: 

- hX^a = -V;r (J//.12.a) 

ya = -hX^i-Vl. (Ill.U.b) 
Using equation (III. 12. b) we eliminate Va in (III. 11) and obtain: 

Ls = l9^JMJ + hU. - V% g,, = hl.hl^, hi = h] + V^bl,, V^^V^- \v^Vl 

(1 1 1. 13) 

The euilibriums of the two systems are the same. In equilibrium Uao — —Vj; (see 
(III. 12. b)), inserting this into equation (III. 12. a) we obtain: 

v,r = - Iv.'v,'),, = vt = o. {in.u) 

Without loss of generality we consider the stability of an equilibrium given by Uao = O,qio, 
where qio are defined by equation (III. 14). 

The energy criterion of the two systems is also the same: For a quadratic system 
such as Lg the energy criterion takes the form Vfj\oSqiSqj > 0, while for a gyroscopic 
system such as Lyj we demand that A^F™ = {dua + V^^^^dqi)^ + {V^-j — V^j^V^.j)dqiSqj = 
(Sva + Va,i3(li)^ + VfjloSqiSqj has a definite sign, i.e. VfjloSqiSqj > 0. 
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It is not to be understood that transforming the weak Lagrangian into Lg has no 
benefit from the point of view of stability analysis. In fact the strong principle provides 
us with necessary and sufficient conditions beyond what we obtain from (III. 8). Every 
perturbation of Ls implicitly includes all dv perturbations, that is we obtain necessary and 
sufficient conditions for more general perturbations than are obtained in the two mode 
analysis of equation (III. 8). For example for a single mode, say 5qi of Lg the gyroscopic 
term vanishes and we have a necessary and sufficient condition of stability: 



this can be compared to a two mode analysis of L^, with 5qi and 5u2 for which we obtain: 



Thus, necessary and sufficient conditions obtained from Lg are more strict. However, in 
Lg one excludes initialy some of the less general perturbations oi 5v in this sense analysis 
through Ls misses some of the necessary and sufficient conditions available in the dynamics. 
We conclude that Lg and Lyj are complementary ways for obtaining necessary and sufficient 
conditions of stability, for sufficient condition the methods give the same results. 

IV. The Energy Criteria for Stability of 2-D Flows 

IV. 1 Basic Identities of Perturbation Theory 

In order to obtain the concrete barotropic flow form of and S'^G of equations 
(III.4.b) that are needed for stability analysis, a few basic identities and notations will be 
given here, this will make our future calculations easier. 



^11 = - KKi > 



(///.15.a) 



(///.15.&) 
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We look at the following Eulerian small displacements: 

aoix.y) ao{x,y)+Sa{x,y), /3o{x,y) /3o{x,y) + 5(3{x,y), Mx^y) iyo{x,y) + 5iy{x,y) 

(IV.l.a) 

and also: 

^co ^co + S^c, bo^bo + 5b. (IV.l.b) 

The subscript o denotes stationary quantites. It is convenient to introduce the often used 
^, which is defined by: 

a{R + i) + 5a{R + i} = a{R) p{R + i) + S^{R + i) = P{R) {IV.2.a) 

to order 1, 

Sa = -i- Va, 5(3 = V/3. (IV.2.b) 
Having defined ^, we can also define the Lagrangian displacement A: 

A = 5 + i-V (/F.3) 

And ofcourse Aa = A/3 = 0. The following identities will serve us at future calculations: 

AR = i, 5V = V5, AV = VA-Vf-V, A^=0. (IVA) 

The variations of both the surface density S and the velocity v can be derived from 
equations (II. 5) and (11.11) (for details see paper I see also paper II): 

AS = -SoV -i, 5E = -V- (EoO {IV.5.a) 

Av = -VC • vq + VAiy, 5v = -Vf • vq + VAiy - ^ ■ Wo, {IV.5.b) 
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we shall also need the second variation of velocity: 

A^v = • Av + VA^v. {IV.b.c) 

IV. 2 The Variation of The Potential V 

We now give the concrete expression of (III. 3) for the potential of a barotropic flow. 
The potential part of the Lagrangian (11.12) is: 

T.(fx-h-{P-PQ)-nc{J -Jq). (/F.6) 

Notice that this is also the energy according to equation (III. 2). 

The condition that AV = for arbitrary perturbations of the form (IV. 1) is that Euler 
and the continuity equations of stationary motion are satisfied in moving coordinates as 
well as the fixation of linear and angular momentum: 

{vo -bo- n^olz X R) ■ Vvo + Qcolz x -Co + V(/io + $o) = 0, (IV.7.a) 

and also: 

V-{J:o{vo-bo-n^olzXR)) = 0, J = Jo, P = Po- {IV.7.b) 

For details see paper I. 

We now look at the second variation 5'^V = A'^V at a configuration satisfying equa- 
tions (IV.7): 

= J {(Aw)2+t/.A2'[;+A[f.V(/i+$)]}|oSod^a;-2A6-AP-2AncAJ-6o-A^P-OeoA2 J, 

(IV.S.a) 
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V 



I 



where: 



and: 



AJ= 1^- J {RxAv+^xvo)^o(fx, AV= 1^- J {2^x Av+RxA^v)i:od'^x. (IV.S.c) 
We demand now that: 

A J = AP = {IV.9) 
this wiU constrain the perturbed configuration to have the same hnear and angular mo- 

— * 

mentum as the stationary one and will define hopefully both A& and Aflc- The perturbed 
potential A'^V will now have the form: 

AV = J{{Avf + v-A'^v + A[i-V{h + ^)]}\oi:od^x-bo-A'^P-ncoA'^J, (IV.IO) 

In configuration having special symmetry QcO or 6o or both, may not be defined by equation 
(IV.T.b), in this case we will demand: A^ J = or A^P = or both. In order to see whether 
the bilinear form A'^V is positive around a certain stationary configuration we will have 
to diagonolize it first this will be done for a few examples in later parts of this paper. 
IV. 3 The Variation of The Kinetic Term G 

Although sufficient conditions of stability can be derived from expression (IV.IO) by 
checking its positivity for certain perturbations, one cannot claim any new knowledge of 
the stability of the system, unless the perturbation is the most general one. The reason for 
this is that if A'^V > it may yet be negative for another form of perturbation and hence 
stability cannot be claimed. If on the other hand A^V^ < the system is not unstable 
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because as we emphasized the /S?V > condition is only sufficient. Thus, in order to 
achive new information we must check how the perturbation under consideration appears 
at the kinetic term 5^G = A^G. As we saw in section III we will have also a necessary- 
condition if the perturbation is of the double mode type. The term G of the Lagrangian 
(11.12) is of the form: 



G = J -a$ - z>] ^(fx 



The first variation of this expression is: 



(IV.ll) 



AG^ 



-aA$ - Av 



{IV.12) 



notice that A/ = (A/) - C • V/ hence: 



AG 



/ 



C-v-{Ai^) 



{IV.13) 



For the second variation near equlibrium we obtain: 



A^G 



/ 



^- A{;+ AC-v- A(Ai/) 



Tiod X. 



(IV.U) 



Using equation (IV. 5. b) we obtain: 



A'G 



2^-Av + {A^u) 



{IV.15) 



The last term is a full time derivative and does not contribute to the equations of motion 



and thus can be neglected, finally we obtain: 



A^G 



-I 



2^-Av 



J^od^x. 



{IV.16) 
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V. Applications to Galactic Flows of Uniform Rotation 

V.l. Stationary Configurations 

The following is taken from Binney and Tremaine (1987) with some changes of nota- 
tions. We consider disks of fluid rotating with uniform angular velocity Q = IgQ and the 
velocity field in inertial coordinates is thus: 

VQ = nxR = VLR^V^p {V.l) 

is the polar angle, R the radial distance. And the density is: 

So = ScSo(i?) {V.2) 

Inserting equations (V.l) and (V.2) into (IV.7) (taking into account the symmetry of the 
potential and the enthalpy derived) we see that the variable VtcO is not defined due to the 
symmetry of the problem and also: 

60 = (K3) 
Pressure and specific enthalpy are given by arbitary equations of state: 

P = P(E), h = h{T). {VA) 

V.2. Global Properties of Perturbed Configurations. 

We now look at the first order perturbation of angular momentum A J of the above 
configurations as it appears in equation (IV. 8. c). We calculate the first term of the above 
equation and use both equations (V.l) and (IV. 5. b): 

j Rx A^loSod^x = j (-Vf • vq + VAz/)Sod^a; = -VL J i- KLQ(fx. {V.b) 
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Calculating the second term we obtain: 
Since (V.b) + {V.6) = 0, we derive that 

A J = {V.7) 

Since Qc is not defined we must demand A^J — 0. Taking J = we are able to obtian 
differential equality convenient for calculating A'^V. It is obtained from (A. 8) in appendix 
A 

A'^V = JiEoiSvf +dJ:d{h + ^)}(fx {V.8) 

V.3. ^ defined by scalar functions 

It will appear very convinient to define ^ in terms of two independant non dimensional 
infinitesimal scalars 77 and ip as follows: 

a is a scale to be defined later. 77 is thus defined in terms of ^ by 

Some boundary conditions are needed to make rj unique. If we take, say, 

v\b = {V.ll) 

equation (V.IO) has a unique solution. Equation (V.IO) represents also the condition of 
integrability of (V.9), considered as a set of two first order differential equations for t/j, 
given ^ and ry. Thus the equations are integrable and define up to a constant. 
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Following equation (IV. 5. b) and equations (V.l) and (V.9), 5v and Aw may be written 
as a gradient plus a rotational, always a convenient form for vector fields: 

5v = a^niVC + rot(2r/l^)] {V.12.a) 



in which 

1 



{Au - ^ ■ vo) - 2i(; {V.13). 



VI. Application to a Uniform Rotating Sheet 

VI. 1. Stationary Configurations 

In this section we study a specific case of the flows described in section V. We consider 
an infinite uniform rotating sheet. Since the density is uniform we have: 

S = 1 (VIA) 

There is no self gravitational field in the plane of the sheet, however, to maintain equilib- 
rium we assume the existance of a radial symmetrical external potential: 

^0 = MR) (^^-2) 
Euler's equations (IV. 7. a) relates to ^o{R): 
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VI. 2. Fourier Decompositions of ^ , 5h, and C, 

In order to diagonolize the second order perturbation of the potential (V.8) around 
our present stationary configuration, we will write our perturbed quantities as fourier 
transforms. First we write r] and t/^ as a fourier transform: 

r]{R) = — / 7^{k)e'^-^d^k {VIA.a) 



27r 



— oo 



oo 



27r 

Next we calculate the Fourier decomposition of Using equations (IV. 5. a) and 
equation (V.9), taking the scale a = 1, (5S can be written as: 

and from equations (VL4) we find that 5T, has the following expansion: 

5T.= — 5^{k)e'''-^d^k {Vl.^.a) 

in which 

5T.{k) = T.ck^r]{k) {Vl.b.h) 
With the expansion of (5E we obtain directly the Fourier decomposition of 5h: 

9h.^ vl 1 r .^.^^„*^.2, 



5h = dJ:{k)e''-''d'k {VI. 6) 

where Vs is the sound velocity. From the varied Poisson equation we calculate the perturbed 
gravitational potential: 

A5^ = ATrGS^ESoiz) {VI. 7) 
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where 5d{z) is the Dirac function, we find the solution which has been given by Binney 
& Tremaine (1987): 



/oo 
SE{k)^e'''-^(fk. {VI. 8) 

-oo 



27rG ^t.s 



The only quantity which we nead to decompose in order to obtain the expansion of Sv is 
C which we now write as: 

C=— C,{k)e'^-^(fk. {VI. 



oo 



27r 

VI. 3. Global Constraints 

According to equation (IV. 9) we must demand that AP = A J = 0. However, A J = 
for all possible perturbations as we derived in equation (V.7). Thus, we must take care 

— * 

only of AP. In addition, due to the symmetry of the problem we saw in section V.2 that 
we must also demand A^ J = . replacing equation (V.12.b) into (IV.S.b), we obtain: 

AP = T.c^ J V(C + ip) + rot{r]l:,)(fx = {VI.IO) 

using stokes theorem this can be written as: 

AP = j>[{C, + -0)1^ xdR + rjl^ ■ dR] = j{C, + -0)1^ xdR = 0. {VI. 11) 

In order to fix the linear momenta we take both tp and ^ to be zero at infinity. 

^ = 0, C = P^oo {VI.12) 

The constraint A^J = becomes, after inserting equations (V.12), (VI. 5) and (IV. 5. c) 
into the second part of equation (IV. 8. c): 



AV = -2Scn / (A-^d^x = {VI.13) 
J dip 
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This constraint can be satisfied by taking r]{k) = r]{k) this will not modify the stability 
condition which which will obtain in the next section since it depends only on k. 
VI. 4- Necessary and Sufficient Conditions for Stability of Infinite Sheets 

A'^V is given by (V.8) near the stationary configuration. To obtain the Fourier decom- 
position of A^F, we replace STj , Sh, 5$ and Sv by their respective expansions (VI. 5. b), 
(VI.6), (VI.8) and (VI.9) with (V.12.b), in A'^V taking account of constraints (V.ll), 
(VI. 12) and integrating we obtain the following results: 

AV = 2Ec J k'^l^'^lCf + (4^^^ + k^v^s - 27rGEc\k\)\rjf]d'^k {VI.U) 

Following section IV.2 the inequality A'^V > 0, calculable from (VI. 14) gives sufficient 
conditions of stability. From (VI. 14) we see that A^ is positive for all k's such that: 

41]2 + k'^v^ - 2TTGEc\k\ > {VI.15.a) 

this means that the sheet is stable to any wave-length if: 

The condition (VI. 15) becomes sufficient and necessary when the gyroscopic term for dy- 
namical perturbations A'^G couples only pairs of modes (see section IV. 3). 

We can make a Fourier decomposition of A^G, using ^ given by (V.9) and and Av 
given by (V.12.b), and using expansions (VI. 4) and (VI.9): 

A^G = 2T,c^ J k^C,'f]*d^k + complex conjugate. {VI. \Q) 

Necessary and sufficient conditions of stability are obtained from (VI. 15) since A^G couples 
only pairs of modes of 77 and C, . The condition (VI. 15) was previously discovered by Binney 

20 



& Tremaine (1987) using the dynamical perturbation method, this was done by our method 
without solving a single differential equation. 

VII. Applications to Maclaurin Disk 

VII.l. Stationary Configurations 

In this section we study a diflFerent case of the flows described in section V. This time 
we study the stability of of the Maclaurin disk which is a finite uniform rotating disk. The 
density of this disk is given by: 



p2 

^ = \I^-^ = X R<Rd, {Vll.l.a) 



D 

S = R>Rd. (Vll.l.b) 

Rd is the radius of the disk. We shall take a = Rd- The self gravitational field in the 
plane of the disk is: 

$0 = -^^0^^ + const., = "^^^ {VII. 2) 
2 2a 

^0 is the angular velocity of a test particle on a circular orbit. See Binney and Tremaine 

(1987) for further details. Here we allow only a certain equation of state such that the 

Pressure and specific enthalpy are given by: 

P = kS^ (yj/.s.a) 

h = ^reS^ {VII.3.b) 
Euler's equations relate Q to Oq, and a: 

Q^ = Ql- {VII. 4) 
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VII. 2. Spherical Harmonic Decompositions of ^ , SH, Sh, and ^ 

We now decompose r) and in normalised spherical harmonic functions of x defined 
in equation (VII.l) and (f. We have: 



oo oo 



^ = E E mmVr{x)e''^'' + c.c. (Vll.b.a) 



l=m m=0 



oo oo 



^ = E E ^imrnx)e''^'' + c.c. {Vll.b.b) 

l=m m=Q 

r]im and i/jim are arbitary, infinitesimal, complex numbers. Associated Legendre polyno- 
mials are defined in the range — 1 < x < 1, but rj and ip are only defined in the range 
< X ^ 1- We may extend rj, t/) to the negative region of x in any way we want. However, 
these functions have bounded gradients at x = 0. Indeed, since 

\dRr]\ < oo, \dRiP\ < oo (VII.G.a) 
1 1 

;d^r]\ < oo, l^d^'iPl < oo {VII. Q.b) 



hold and since 



9ii = -— yr^a^, {VII.l) 

ax 



then. 



5x^lx=o = ^xV'lx-o = 0. {VII.%) 
Equations (VII. 8) will be satisfied if we make symmetrical continuous extensions 

^(X) = ^(-X), {VII.9.a) 

V'(x) = V^(-x), {Vll.^.h) 
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whose expansions are given by (VII. 5) with (l — m) even. The expansions in the domain 
< X ^ 1) of arbitrary continuous rj and ip that satisfy (VII. 6) everywhere are then given 
by (VII. 5) with (/ — m) even [Arfken (1985)] . The boundary conditions (V.ll) gives the 
foUowing relations among the r/^m's: for every m > 0: 

oo 

J^VimVPiO) = 0, l-m even (F/LIO) 

l=m 

Equations (VII. 10) define, say, rju in terms of rji^mj^l- 

We can now calculate the spherical harmonic decomposition of ST, using equations 
(IV.5.a) and (V.9): 

SE = -a^Ec[V ■ (xV?7) + Vx x V^]. {VII. 11) 

and from equation (VII. 11) and (VII. 5) we find that ST, has the following expansion: 



oo oo 



5T = TcJ2J2^i^ ^'""^ + I -m even {VII.12.a) 

l=m m=0 ^ 



in which 



S/m = + 1) - ■m^]ilim + imtpim = kimVlm + imtpim {Vll.U.b) 

With the expansion of (5E we obtain directly through (VII. 3. b) the spherical harmonic 
decomposition of Sh: 



5h = 3kTI J2 Y1 ^irnVrixV'^^ + c.c. I - m even {VII.IZ) 

l=m m=0 

and from the varied Poisson equation : 

AS^ = AttGSTSd (z) {VII. 14) 
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where 5d{z) is the Dirac function, we find the solution 6^ which has been given by Hunter 
(1963): 



oo oo 

5^ = rtla^ J2 Yl ^«mPr(x)e*""^ + c.c. I - m even (VII.15.a) 

l=m m=0 

in which 
where 

9im= y^.TJ'^!?,"^^:.^ <1 ^>0, I -m even {VII.15.c) 
We now expand Q in spherical harmonics: 

oo oo 

C = E E Om7'r(x)e'"^^ + C.C., I - m even 

l=m m=0 

VII. 3. Global Constraints 

— * „ 

In this section as in VI. 3 we must take care of AP = and A J = . Inserting 
equation (V.12.b) into the first expression of (IV. 8. b) and using expansions (VII. 5) and 
(VII. 16) we obtain the result: 

Cii = -V'ii-%i- {VII.17) 

i.e. in order to fix the momentum we must take (n to be a dependent variable. The 
constraint J = becomes, (the calculation is in appendix B): 

A2j = 47rSca4l7 2^m[iO^S*i„+m(|V'/mP - \mm\^)] + c.c. = {VII. 18) 



{VII.16) 
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VII. 4- Necessary and Sufficient Conditions for Stability of Maclaurin Disks 

In order to study the stability features of the Maclaurin Disks we shall have to di- 
agonolize A'^V given by (V.8). To obtain the spherical harmonic decomposition of A^V, 
we replace 6T, , 6h, by their respective expansions (VII. 12), (VII. 13) and (VII. 15). Also 
we write Sv by the represantation (V.12.a) using the expansions (VII. 5. a), (VII. 16) . The 
result is as the follows [detailed calculations are given in appendix C] : 

A V = 47rScJ^^a^ ^hm - ^)\riim? + kim\Clm + T—rilm?+ 

[-1 + (1-^,^)^]|S^^|2 (1/J/.19) 

where we have taken into account the constraint (VII. 17). Following section IV. 2 the 
inequality A'^V > 0, calculable from (VI. 19) with the constraint (VII. 18), gives sufficient 
conditions of stability. A'^V has the form fpA'^ — B > 0. If S < 0, the disk is stable for 

q2 

any ttj- For instance, Maclaurin disks are stable to any perturbation that keeps the same 
densities at the displaced points {ST, = 0). We are naturally interested in perturbations 
that might upset stability and for which B > 0. If 5 > then must satisfy the following 
inequality 



Q=7v?< 



^0 S(ILm)>l,m=0 l^imP - ^{hm - ■^)|^«mP - hm\Clm + |^?7«mP 

(V^//.20) 

The condition (VII. 20) supplied by the constraint (VII. 18) becomes sufficient and necessary 
when the gyroscopic term for dynamical perturbations A^G couple only two modes (see 
section III). It is therefore important to calculate A^G. We can make a spherical harmonic 
decomposition of A'^G given in equation (IV. 16). A straightforward subtitution of ^ given 
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by equations (V.9) and (VII.5) and given by equations (V.12.b), (VII.5) and (VII.16) 
leads to the following expression for A^G: 

A^G = 47rScJ^a^ [^imCim + i'l^iVimVim - i^im'^lm)] + complex conjugate {VI 1. 21) 

l,m^l,l 

in which (VII. 17) and (VII. 18) has to be taken into account. Necessary and sufficient 
conditions of stability are obtained from (VII. 20) when A^G couples pairs only. 

VII.5. Symmetric and Antisymmetric Single-Mode Perturbations 

In a single-mode analysis condition (VII. 18) reduces to: 

m[iC/mS *im +m{\^im\^ - \r}im\^)] + c.c. = {VII. 22) 

For symmetrical and antisymmetrical modes, we take either real or imaginary components 
of riim and i'^im in the spherical harmonic expansion, so the Fourier expansion contains 
either cos{nnp) or sin{m(^). For such perturbations, A^G becomes: 

A^G = 47rSc^^a^ ^ [^imCim] + complex conjugate {VII. 23) 

where we have taken into account the constraint (VII. 17). Notice that we have used without 
loss of generality we assume that Him is real (and so is rjim and iipim)i and equation (VII. 23) 
becomes: 

A^G = 87rScl^a^ J2 PimOmi?]- {VI 1. 24) 

We now have a pair coupling of Tjim and (imR, hence we can achive a necessary and 
sufficient condition of stability. The imiginary part of C«m can now be obtained through 
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equation (VII. 22), assuming non-radial modes (for radial modes equation (VII. 22) is triv- 
ially satisfied): 

Qmi = m =^ {Vll.Zb) 

It is advantageous to introduce the following notations: 

z = — fc/m-^r~" — oo < 2; < oo {VII.26.a) 

and 

a:=T2- 0<x<l. {VII.26.b) 
For one single mode (Z, m), equation (VII. 20) with equation (VII. 25) reduces to in term of 



Q<Qlm= ^V/t^ {VII.21.a) 

we have assumed that ClmR is equal to zero since any other value of ClmR will make Qim 
bigger, and we are intrested only in its lowest value. P4{z) is a polynomial of order 4 in z, 

P^{z) = 1 - T^{(- - 1)[(1 - x)z'^ + 4(1 - x)z^ + Qz^ + Az] + -} {VII.27.b) 

Pi{z) has one and only one real maximum for any x. We are only interested in values of z 
for which ^4(2;) > 0. For ^4(2;) < 0, A^F > for any value of Q. The maximum of Pa^z) 
is obtained for: 

i_ 

z^ax{x) = -^L=[{^/T^ - 1)3 + (v^r^ + 1)^] - 1 (yi/.28.a) 

— X 

for which 

P4(^)ma. = l-T:^{(i-l)[(l-a;)4„, + 4(l-a;)2;i,, + 64,, + 4w] + ^} = 1-^ 

Kim X X l^lm 

{VII.28.b) 
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To any pair of values {l,m) corresponds a value x defined by (VII. 26. b) and a value y{x). 
Following (VII. 27. a) and (VII. 28), we must thus have 

Q < {Qlm)min = ^^^^ {Vn.29) 

The smallest minimum of Qim is obtained for (Z,m) = (2,2) for which {Q22)min = |- 
Therefore the necessary and sufficient condition of stability with respect to symmetric or 
antisymmetric single mode perturbations of Maclaurin disks is 

Q<\ (y//.30) 

VII. 6. Comments 

a) Binney and Tremaine (1987) have given the following dispersion relation for the uj - 
modes of dynamical perturbations: 

cu^ - cUri^Q'^ + kimi^lil - Qim) - J^^]} + 2mVL[VLl{l - giro) - f^^] = 0, Ur = U- mfl 

(VII.Sl) 

Stability holds if the non spurious co^s are real roots. Equation (VII. 31) is a third order 
polynome. The condition for a polynome of the form + a2X + as = to have only real 
roots is given in standard handbooks [for instance: Schaum's Mathematical Handbook 
(1968)] 

(|)3+(|)2<0 {VII.32) 

Inequality (VII.32) for equation (VII. 31) is exactly our inequality (VII. 29). The results 
of Binney & Tremaine were obtained by solving three couple differential equations, our 
results were obtained without solving a single equation. 
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b) For radial modes, m = 0, , J = and 

A^G = 47rScna4^[E;oOo]. {VI 1. 33) 

1=0 

that is all radial modes are only pair coupled. Moreover, since [see equation (VII. 12. b)] 

S^o = /(/ + l)mo (F//.34) 

all Z-modes are decoupled, A^V is a sum of squares of independent modes. In this case 
one obtains stability limits for non single modes. The most unfavorable limit, however, is 
again that given by Qio for which y = 4 and ku = 1(1 + 1): 

Q<Qlo = T^^^ a>2), 910= ^,J'., {vn.35) 

c) Notice the similarity of the mathematical analysis in the weak energy method and the 
strong method of paper II. 

VIII. Application to a Rayleigh's Shear Flow 

VIII. 1. Stationary Configurations 

The problem of the stability of a shear flow was considered by most developers of 
energy methods such as Arnold (1969), Holm at al. (1983) and Grinfeld (1984). It inter- 
esting therefore to compare it to our own weak energy method. We consider a flow in a 
rectengular pipe, which occupies the domain: [0 — L, — 1] that is the rigid walls of the 
pipe are the lines j/ = 0, y = 1 while the lines a; = 0, x = L are identified, in this way we 
construct a "torus". x,y are the Caretesian coordinates. We define the fiow such that: 

vo = W{y)l,. {VIII. I) 
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The uniform density is set to unity 

So = 1 (yill.2) 
Pressure and specific enthalpy are given by arbitary equations of state 

P = P(E) (Vlll.S.a) 

h = h(T,) (Vlll.S.b) 

There is no potential of any kind in this problem. Inserting equations (VIII. 1-3) into the 
equations of motion (IV. 7) we obtain the following: 

= 0, boy = (VIIIA) 

box is left undefined as was flc in the models considered before. Notice in this problem we 
do not have a "free boundary" but a fixed one, this modifies the possible perturbations 
one can construct such that: 

Su{0, y) = Su{L, y), ({0, y) = i{L, y) {Vlll.h.a) 

U^,Q)^iy{x,l)^Q {Vlll.h.b) 

since we consider here a torus with rigid boundaries. 

VIII. 2. Global Constraints 

In this case fic = 0, b^y — Q are defined but not feoa;- This implies that not only do we 
have to fix: AP = 0, A J = but A^P^; = as weU. Inserting equations (VIII.l), (VIII.2) 
and (IV. 5. b) in the first equation of (IV. 8. b) we obtain: 

AP = y Av(fx = j {-Wi^W + WAu)(fx = ly j -^W + ^^d'^x = (VIII.Q.a) 
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Notice that the x-component which contained only boundary terms vanished do to equa- 
tions (VIII. 5). The second order A^P^; = becomes using equation (IV. 5. c) 



A^P, = J A\a:d^x = -2 J -^■Av + —Q^d^x ^ -2 J ■ Avd^x (VIII.G.b) 
A^u term vanishes due to integration by parts. The angular momentum is fixed by: 

A J = 1^ ■ y ixvo+Rx Avd^x = J -^yW + l^-{Rx Av)(fx {VIII.Q.c) 

VIII. 3. AV 

We start from equation (IV. 8. a). Using equations (A. 6), (VIII. 3) we have: 

^ ^ ^ ^ dh ^ ^ 

A(e • V/i) = e • V(5/i = e • V(^|o5S) = i ■ V{vl5T) {VIII. 7) 

we also used the definition of the velocity of sound Vg. This means that the "potential" 
part of equation (IV.S.a) becomes using equations (VIII. 5), (IV. 5. a) and integrating by 
parts: 

A% = j A[i-Vh]d'^x^ Jv^JE'^d^x {VIII.7) 

Using equation (IV.3) and (VIII. 1) we have: 

Av = Sv + ^yW'la;. (VIII. 8) 

Putting equations (VIII. 8), (IV.5.c) into the "kinetic" part of equation (IV.8.a) we have: 

AVfc = J [Sv'' + 2^yW'Sva, + ^IW" - 2W^ ■ {5v + ^yW%)]d^x {VIII. 8) 

A?v term vanishes due to integration by parts. We use equation (IV.S.a) to subtitute 
= — 5S — in equation (VIII. 8) and obtain: 



AVfe = y \8v'' + 2iyW'dv^ + ilW" + 2W{W'iy + 8v^){8Y. + ^) - 2W^^H]^'a; 

(F///.9) 
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Integrating by parts the second and last terms and collecting some terms together we arrive 
at the following form: 



The term 6uj = (^^^ — ^|^) which is the variation of the vorticity can be written using 
equations (IV. 5) and (VIII.l) as: 

Soj = ^yW" - SEW' (VIII.ll) 

Inserting equation (VIII. 11) into (VIII. 10) and canceling the apropriate terms we arrive 
at the simpler expression: 

AVfc = J [Sv^ + {5va; + W5Ef - W^^^E^ + WW"^l]d^x. {VIII.12) 

Thus: 

= AVfe + AVp = J [Svl + {5v^ + WSJ:f + {vI-W^)SJ:'^ + WW"^l]d^x. {VIII. 13) 

Prom this we extract the following famous sufficient condition: 

> and WW" > 0. {VIII.U) 

This was obtained by Lord Rayleigh (1880), Arnold (1966), Holm et al. (1983) and Grinfeld 
(1984). This is obiously a weak condition since we did not take full advantage of the 
constraints (VIII. 6). 

The author would like to thank J. Katz for critical discussions and for his suggestion 
to add section (III.l). 
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Appendix A: A'^V Variational Identity (V.8) 

We start from (IV. 10). Inserting equation (V.l) into equation (IV.IO) we obtain: 

AV = J{{Avf + {nxR)-A^v + A[i-V{h + ^)]}\o^od^x (A.l) 

Since we demand that A^J — 0, we see from the second part of (IV. 8. c) that we can get read 
of the A'^v term thus obtaining: 

= J {(Avf - n ■ (2|*x Av) + A[i- V{h + $)]}|oSorf^a; (A.2) 

The terms containing Av can be written as a difference of squares: 

(Avf-il-{2ixAv) = (Av-ilxiy-(nxi)^ = (Av-i-Vvof-n^i" = (Syf-n^i" (A.3) 

In the last equahty we have used the second part of equation (IV. 5. b). So equation (A. 2) with 
(A. 3) can now be written as: 

= J {{5vf - + A[f V(/i + $)]}|oSod'x (^.4) 

Using the third operator identity of equation (IV. 4), we may rewrite the h + ^ term in (A.l) 
as follows: 

A[f V(/i + $)] =fVA(/i + $)-f-Vf V(/t + $) =fW(/i + $) + fVV(/t + $)-f (A.6) 
But from Euler's equations (IV. 7. a) we see that 

V{h + ^)\o = -vo-Vvo {A.7) 

Subtituting (A.7) into (A. 6) and A[f • V(/i + $)] back in (A.4), we find that 17^^ cancels out. 
Finally with ST, given in (IV. 5. a), and with some integration by parts, we obtain the following 
much simpler form for A'^V which is that written in (V.8): 

AV = J{J:{6vf + dJ:d{h + ^)}\o(fx (A.8) 
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Appendix B: Spherical Harmonic Decomposition of J 

Following the second part of (IV. 8. c): 

AV =1^. J{2(x Av + R X AV)|oSod^a;. {B.l) 

Inserting equation (IV. 5. c) into (B.l) we obtain: 

AV = 1^- J{2ix Av-2Rx Vf • Av)\oJ:od'^x. {B.2) 

Inserting (V.12.b) and (V.l) into (B.2) gives 

AV = 2a^m^ ■ J{[^V + rotjf] x [V(C + ^) + rot{r]l^)] 
-Rx V[V77 + rotif] ■ [V(C + -0) + rot{ril^)]}i:o(fx. (B.S) 

The spherical harmonic decomposition of A^J is obtained by replacing in equation (B.S) C by 
(VII.16), So by Sex (equation (VII.l)) and r/, ip by (VII.5). Using also (VII.12.b) for S/^, 
the result comes out as follows: 

AV = 47rSca^l^ ^m[iOmS *im +m{\ij;im\^ - \mm\'')] + c.c. (BA) 

and demanding A^ J = we obtain equation (VII. 17). 
Appendix C: Spherical Harmonic Decomposition of A'^E 
Let us start from (A. 8) and take A^J — 0. We shall set 

AV = AVfe + AVp (C.l.a) 



in which 



A^Vk = J[E{Svf]\od''x {C.l.b) 
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and 



AVp = J [dJ:{6h + (5$)] \o(fx. (C.l.c) 
Inserting (V.12.a) in (C.l.b) gives: 

AVfc = a'^n^ J {(VC)^ + 2V • [Crot{2r]l,)] + [rot{2r]l,)]'^}i:od^x (C.2) 

Using Gauss's theorem, we can integrate (C.2) by part, and with SqIb = 0, obtain: 

Vfc = a^n^ J^~^^ ■ (^oVC) - 2CVSo • rot{2r]l^) - 4r7V • (EoVr/)]d\. (C.3) 



^2. 



The spherical harmonic decomposition of A^14 is obtained by replacing in equation (C.3), C 
by (VII. 16), So by T,cX (equation (VII. 1)) and ry by (VII. 5. a). The result comes out as follows: 



AVfe=47rScl^V Yl mim-^)\Vln.\' + klm\Cln. + ^mm\'] (CA) 
l=m,m=Q 

Notice that these are already some of the terms of A^F given in (VII. 18) 

The spherical harmonic decomposition of /Sp'Vp follows by inserting in (C.l.c) the respective 
expansions of 5S in (VII.12), 5h in (VII.13) and 5^ in (VII.15). For A^V^p, we obtain : 



AVp = 2TTa^nl E + 3«:ScSr^] + c.c. (C.5) 

l=m m=0 

With (Vll.lS.b) and (VII.4), A^F^ can also be written: 

A = 47rEca^l72 J] [_1 + (1 _ ^^^)_0]|e,^|2 (c.g) 
The sum of A^V^ given by (C.4) and A^ii^p of (C.6) is the expression of A^v written in (VII. 18). 
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